COSMOS* morphological classification with ZEST 
(the Zurich Estimator of Structural Types) 
and the evolution since z = 1 of the Luminosity Function of early-, 

disk-, and irregular galaxies 

C. ScarlataS C. M. CarolloS S.J. Lilly^, M. T. SargentS R. Feldmann^ P. Kampczyk^, C. 

Porciani^, A. Koekemoer^, N. Scoville^, J-P. Kneib^, A. Leauthaud^, R. Massey^, J. 
Rhodes2'^ L. Tasca^, P. Capak^ C. Maier\ H. J. McCracken^ B. Mobasher^, A. Renzini^ 
Y. Taniguchi^ D. Thompson^'^, K. Sheth^'^o^ M. Ajiki", H. Aussel^^.is^ rp Murayamai\ D. 
B. Sanders^^^ S. Sasaki^l•l^ Y. Shioyal^ M. Takahashi^^-i^ 



^Department of Physics, Swiss Federal Institute of Technology (ETH-Zurich), CH-8093 Zurich, Switzer- 
land 

^Space Telescope Science Institute, 3700 SanMartin Drive, Baltimore, MD 21218 
^CaHfornia Institute of Technology, MC 105-24, 1200 East CaUfornia Boulevard, Pasadena, CA 91125 
''Laboratoire d'Astrophysique de Marseille, BP 8, Traverse du Siphon, 13376 Marseille Cedex 12, France 
^Jet Propulsion Laboratory, Pasadena, CA 91109 

^Institut d'Astrophysique de Paris, UMR 7095, 98 bis Boulevard Arago, 75014 Paris, France 

^Dipartimento di Astronomia, Universita di Padova, vicolo dell'Osservatorio 2, 1-35122 Padua, Italy 

^Subaru Telescope, National Astronomical Observatory of Japan, 650 North Aohoku Place, Hilo, HI 
96720 

^Large Binocular Telescope Observatory, University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721 

^^Spitzer Science Center, California Institute of Technology, Pasadena, CA 91125, USA 

^^Astronomical Institute, Graduate School of Science, Tohoku University, Aramaki, Aoba, Sendai 980- 
8578, Japan 

^^Institute for Astronomy, 2680 Woodlawn Dr., University of Hawaii, Honolulu, Hawaii, 96822 
^■^Service d'Astrophysique, CEA/Saclay, 91191 Gif-sur-Yvette, France 

^"^Physics Department, Graduate School of Science and Engineering, Ehime University, Japan 

* Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA Inc, under NASA contract NAS 5-26555; also based on 
data collected at : the Subaru Telescope, which is operated by the National Astronomical Observatory of 
Japan; Kitt Peak National Observatory, Cerro Tololo Inter- American Observatory, and the National Optical 
Astronomy Observatory, which are operated by the Association of Universities for Research in Astronomy, 
Inc. (AURA) under cooperative agreement with the National Science Foundation; and the Canada-France- 
Hawaii Telescope operated by the National Research Council of Canada, the Centre National de la Recherche 
Scientifique de France and the University of Hawaii. 



-2- 



ABSTRACT 



Motivated by the desire to reliably and automatically classify the structural 
properties of hundred-thousands of COSMOS galaxies, we present ZEST, the 
Zurich Estimator of Structural Types. 

To classify galaxy types and structure, ZEST uses: (i) Five non-parametric di- 
agnostics, i.e., asymmetry A, concentration C, Gini coefficient G, 2nd-order mo- 
ment of the brightest 20% of galaxy pixels M20, and ellipticity e; and {ii) The 
exponent n of single-Sersic fits to the two-dimensional surface brightness distri- 
butions. To fully exploit the wealth of information while reducing the redundancy 
present in these diagnostics, ZEST performs a principal component {PC) Anal- 
ysis. We use a sample of ~ 56, 000 Iab < 24 COSMOS galaxies to show that 
the main three PCs fully describe, without significant loss of information, the 
key aspects of the galaxy structure, i.e., to calibrate a three-dimensional ZEST 
classification grid of axis PCi, PC2, a nd PC^. We dern onstrate the robustness of 
the ZEST grid on the z = sample of iFrei et al.l (119961 ). The ZEST classification 
breaks most of the degeneracy between different galaxy populations that affects 
morphological classifications based on only some of the basic diagnostics included 
in the ZEST structural analysis. 

As a first application of the ZEST structural classification scheme, we present 
the evolution in the 0.2 < z < 1 redshift regime of the Luminosity Functions 
of COSMOS galaxies of different structural types (i.e., early, disk and irregular 
galaxies) and, for disk galaxies, of different bulge-to-disk ratios. Overall, we find 
that the bright-end of the luminosity function up to a redshift 2 = 1 is consistent 
with a pure-luminosity evolution (of about 0.95 magnitudes at 2; ~ 0.7). We 
highlight however two trends, relative to the local Universe, that are in general 
agreement with a down-sizing scenario for galaxy formation, i.e.: (1.) A deficit of 
a factor of about two at 2; ~ 0.7 of Mb > —20.5 structurally-classified early-type 
galaxies; and (2.) An excess of a factor of about three, at a similar redshift, of 
irregular galaxies. 



Subject headings: cosmology: observations — cosmology: large scale structure of 
universe — cosmology: dark matter — galaxies: formation — galaxies: evolution 
— surveys 
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Introduction 



Galaxy morphologies are strongly correlated with the star formation histories of galaxies, 
are a proxy for their dynamical structure, and are thus a key diagnostic of galaxy evolution. 
As state-of-the-art surveys of faint galaxies at (low and) high redshifts such as COSMOS 
now return millions of galaxies, it becomes essential to have robust tools for automatic 
morphological and structural classifications. 

Motivated by the need to automatically derive a quantitative description for the mor- 
phology of the large number of COSMOS galaxies, we present the Zurich Estimator of 
Structural Type (ZEST), a powerful classification scheme that combines: 
— A Principal Component Analysis (PCA) of five non-parametric diagnostics of galaxy 
structure, i.e., asymmetry A, concentratio n C, Gini coefficient G, 2nd-order moment of 
the b rightest 20% of galaxy pixels M20 (e.g 



Abraham et al. 2003: Lotz. Primack &: Madau 



2004J ). and the ellipticity e. The PCA reduces the redundancy of information that is present 
in these five diagnostics, without substantial loss of information on the galaxy structure. 
— A parametric description of the galaxy light, i.e., th e exponent n of a single-Sersic fit to 
its surface brightness distribution (ISargent et al. II2006I . this volume). 

Specifically, ZEST classifies galaxies on the basis of their location in the three-dimensional 
space defined by the main three eigenvectors PCi, PC2 and PC3 that contain most of the 
variance of the original non-parametric quantities, while the n- Sersic indices are used to 
refine the classification of disk galaxies by splitting these in four separate bins of a bulgeness 
parameter (that is correlated with the bulge-to-disk ratio). 

We calibrate the ZEST grid with a sample of ~ 56000 Ub < 24 COSMOS galaxies se- 



duced from the Cycle 12 HST ACS F814W images fIScoville et al. 



(Scoville et al. 


2006; 


Leauthaud et al. 


z = sample of 


Frei et al.l ( 


1996). The 



ZEST grid assigns, to each galaxy: 

(1) A morphological type (=1, 2 or 3 for early-type, disk and irregular galaxies, respectively); 

(2) A bulgeness parameter for the type 2 disk galaxies, which coarsely correlates with the 
bulge-to-disk ratio (B/D); the disk galaxies are split in four bins, from pure disk galaxies 
(T = 2.3) to bulge dominated disks (T = 2.0); 

(3) An elongation parameter (four bins, respectively from = for face-on to = 3 for edge-on 
galaxies) ; 

(4) An irregularity and a dumpiness parameter for types 1 and 2 galaxies, which indicate 
whether the galaxy light distribution is regular/smooth or distorted/clumpy, respectively. 

Although the size of the galaxy does not enter d irectly in the ZEST classification, a 
measurement of the Petrosian radius (jPetrosian Ill976l ) is produced as a by-product of the 
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ZEST classification procedure, and is thus also available for all the ZEST-classified galaxies. 



The calibrated ZEST grid in the PC1-PC2-PC3 space can be used to automatically 
classify galaxies according to their structural properties. 

While different combinations of some of the basic diagnostics that we use to construct 
the ZEST grid have been extensively used in the literature to classify galaxy types, these 
combinations are effectively "projections" on lower-dimensionality planes of the ZEST PCi- 
PC2-PC3 space. These projections mix very different galaxy populations (e.g., elliptical 
galaxies with regular edge-on disks in the M20-G and C-M20 planes, etc.), and thus lead to 
heavy contamination of galaxy samples and erroneous statistical assessments concerning the 
evolution of specific galaxy populations. The three-dimensional PC1-PC2-PC3 space breaks 
the degeneracy between these different galaxy populations. 

As a first application of the ZEST structural classification of COSMOS galaxies, we 
present the redshift evolution since 2; = 1 of the Luminosity Functions (LFs) of galaxies with 
different structural types (including different "bulgeness" for the disk galaxies). In order 
to anchor the COSMOS measure ments to the loca l Universe, we use the sample of 1813 



Sloan Digital Sky Survey (SDSS; lYork et al.l 120001 ) galaxies constructed and discussed by 



Kampczyk et al. (2006). A detailed discussion of the LFs of more elaborated selections of 
gala xy samples that are st ructurally classified with ZEST will be presented elsewhere (see. 



e.g.. iScarlata et al. Il2006bl . for an extensive discussion of the evolution of the LF of elliptical 
galaxy progenitors). 

This paper is structured as follows. In Section 2 we describe the methodology at the basis 
of the ZEST classification: in particular. Section 2.1 describes the non-parametric quantifi- 
cation of galaxy structure. Section 2.2 the principal component analysis in the A, C, G, M20, 
and e space. Section 2.3 presents the calibration of the ZEST morphological grid, and the 
use of the n-Sersic index values to establish the B/D ratio of disk galaxies (Section 2.3.2). 
The performance and reliability of ZEST are discussed in Section 3, where we present the 
test performed on a z = galaxy sample (Section 3.1), and a discussion of the error bars in 
the ZEST structural classification (Section 3.2). In Section 4 we show the advantage of the 
ZEST classification scheme over other approaches. These comparisons show that the ZEST 
classification is substantially more powerful, relative to other schemes that are based on only 
a few of the basic diagnostics that are contained in the ZEST grid, to separate galaxies 
with different structural properties. In Section 5 we present the Luminosity Functions of the 
ZEST-classified COSMOS galaxies, the corrections applied to the LFs, and the sources of 
error. The paper is briefly summarized in Section 6. 
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Throughout the paper we use the cosmology that has been adopted throughout this 
volume, i.e., = 0.25, flm + flA = 1 , and Hq = 70 km s~^ Mpc~^. All magnitudes are 
expressed in the AB system (lOke Ill974j ). 



2. Methodology 

ZEST is based on a combination of non-parametric and parametric quantification of 
galaxy structure, and a PCA to reduce the number of variables. We describe below the steps 
that we followed to construct the ZEST classification grid. 



2.1. Non-parametric quantification of structure 



The ZEST PCA uses, as basic quantities, four widely-adopted non-parametric measure- 
ment s of the light distribution in galaxies (e.g., [Abraham et al.ll2003l : lLotz. Primack fc Madau 
20041 ). plus a measurement of the ellipticity of the galaxy light distribution. Specifically, 
ZEST uses: 

1. ) The concentration C (= 51og^; with rgo and the radii including 80% and 20% of 
the total galaxy light, respectively), which quantifies the central density of the galaxy light 
distribution; 

2. ) The asymmetry A, which quantifies the degree of rotational symmetry of the light 
distribution. A is measured by calculating the normalized difference between the galaxy 
im age and the image rotated by 180°. A correction for b ackground noise is also applied (as 



m 



Consehce et al.ll2000l : iLotz. Primack fc Madau II2004I ) 



I.e.: 



A 



E 



x,y \h^,y) ~ -^180(x,y) 



-B 



180 1 



2 E \h=^,y)\ 

where / is the galaxy fiux in pixel (x, y), /i8o is the image rotated by 180° about the galaxy's 
central pixel, and Sigo is the average asymmetry of the background; 

3.) The Gini coefficient G which describes how uniformly the fiux is distributed among 
galaxy pixels. The Gini statistic assumes values from (if the galaxy light is homogeneously 
distributed among galaxy pixels), up to 1 (if all the light is concentrated in one pixel, 
regardless of its position in the galaxy). Specifically, after ordering the pixels by increasing 
fiux value, G is given by: 
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1 " 



wher e n is the number of pixels assigned to a galaxy, and X is the mean pixel value (iGlasser 



19621); 



4-.) M20 (= logX] Mi/Mtot, with ^ fi < 20% and Mtot the total second order moment), i.e., 
the second order moment of the brightest 20% of the galaxy flux. For centrally concentrated 
objects, M20 correlates with the concentration C; however, M20 is also sensitive to bright 
off-centered knots of light; 

5.) T he ellipticity e = 1 — b/a of the light distribution, as measured by SExtractor (version 



2.4.3. iBertin fc Arnoutslll996l ). SExtractor computes the semi-major axis a and semi-minor 
axis b from the second order moments of the galaxy light; specifically, a and b are the 
maximum and minimum spatial variance (rms) of the object, along the direction 6 where 
the variance is maximized: 



a^ = ^ + \/^^^^^^ + -y, (3) 



6^ = :^- + (4) 



where the second order moments x2, y2, and xy are given by: 



X' = - x', (5) 



(6) 



Z^j liXiyi — 
xy = xy. (7) 

The values of a and b are thus representative of the galaxy ellipticity at large radii, and 
are rather insensitive to details in the internal structure (e.g., bars in disk galaxies, knots of 
star formation, etc). 
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2.2. Principal Component Analysis 

The quantities above provide complementary, but also redundant information on galaxy 
structure. We therefore performed a PCA using the measurements of A, C, M20, G and 
ellipticity e as basic variables. 

The PCA is a classical statistical method for multivariate analysis, which reduces the 
dimensionality of a dataset without a significant loss of information. This is done by trans- 
forming the observed variables into a new set of orthogonal variables, the principal compo- 
nents ^^PC", (with i = 1, ..,n, and n the number of basic parameters, i.e., variables). The 
PCi are ordered so that the first few of them retain most of the variance present in the 
original dataset. The principal components are a linear combination of the original vari- 
ables, and define a new coordinate system obtained by rigid rotation of the original space. 
In the new system, the axes represent the directions of maximum variability in the original 
n— dimensional distribution of points. 

In detail, the dataset is described by an n x m data matrix, (n = 5 in the current 
version of ZEST), and m is the number of galaxies with measured basic parameters. All 
variables are standardized before performing the analysis by subtracting their median value 
(indicated with the subscript M in the following expressions) and normalizing them with their 
standard deviation. Therefore the five variables considered in the ZEST PCA are defined as 
xi = {C - CM)/crc, X2 = (M20 - M2o,m)/o-m2o, X3 = {G - Gm)/ctg, X4: = {e- eM)/cr^, and 
x^ = {A- Am) jo A- 

The directions of the principal components are derived by calculating the eigenvectors 
of the X n covariance matrix of the Xj variables [Sij = ((xj — {xi)){xj — (xj)))). The matrix 
S" > is real and symmetric. Thus, it admits real, positive eigenvalues Aj. By sorting the 
eigenvectors in order of decreasing values of the eigenvalues, an ordered orthogonal basis is 
obtained, with eigenvectors aligned along directions of decreasing variance (A,/ J2j -^j) ^^e 
data. The first few principal components that account for most of the power, i.e., most of 
the total variance Aj- in the dataset, are then used to replace the original n variables 
without any significant loss of information. 



2.3. ZEST calibration with 56000 Iab < 24 COSMOS galaxies 



We calibrate the ZEST classification grid on a sample of ~ 56000, Iab < 24, COSMOS 
galaxies detected in the 260 A CS F814W images acquired during the HST Cycle 12 observing 
period (IScoville et al. II2OO6I ). The total area covered by this fraction of COSMOS is 0.74 



deg2. Details on the COSMOS sample are given in Appendix lA.ll 
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For each COSMOS galaxy, we measured the basic non-parametric quantities described 
in Section 12.11 by computing them on the galaxy pixels (defined using Petrosian apertures, 
see Appendix IA.2p . Figured] shows the behavior of each basic non-parametric diagnostic as a 
function of the others. The contours in each panel enclose respectively 30%, 80% and 98% of 
the COSMOS galaxies in our sample. Global correlations are known to exist between various 
non-parametric coefficients. For example, relatively tight correlations exist between G, M20, 
and C, with objects with high C tending to have low M20 and high G. Any value of C is 
observed for small values of A, while high values of A are preferentially observed in low-C 
galaxies. These trends have already been noted in the literature, and indeed highlight the 
redundancy of information present in these diagnostics. As expected, the ellipticity e does 
not correlate with any of the other parameters, except for a mild positive correlation with 
the concentration for e > 0.6. This is a geometric effect, since edge-on galaxies preferentially 
have high C values. 

The results of the PCA on the normalized COSMOS dataset are presented in Table [1] 
In particular. Columns 2-6 refer to the five principal components derived in the analysis. 
The first row gives the eigenvalue (i.e., variance) of the data along the direction of the 
corresponding PC. The second row shows the fraction of the variance that is explained 
by each of the PCs, i.e., the fraction of the "power" that is contained in each PC; the 
third row lists the cumulative fraction of the variance. In the bottom part of the Table, 
each column lists the weights assigned to each input variable (listed in Column 1), in the 
linear combination that gives the direction of the specific principal component (e.g., PCi = 
-0.54 xxi + 0.60 XX2 - 0.56 x 0:3 + 0.20 x X4 + 0.02 x 0:5). 

In Figure [2] we show the fraction of variance as a function of the corresponding principal 
component. Solid circles refer to the PCA applied to the COSMOS galaxies in the considered 
Iab < 24 sample. Open squares refer to the same analysis, performed however on only the 
COSMOS galaxies brighter than I^b = 22.5. The horizontal line represents the value of the 
eigenvalues that would be expected if the five variables were uncorrelated. The comparison 
of the results for the Iab < 22.5 and the Iab < 24 samples demonstrates the stability of the 
analysis down to the faintest magnitudes in the sample. 

As discussed in Appendix [Bl the scatter in the measured variables that define the PCs 
increases with increasing magnitude (the scatter at Iab = 24 is twice the scatter measured 
at Iab = 22). While the larger scatter in the data could potentially wash out correlations 
present in the original parameters, this effect is negligible down to the considered Iab = 24 
magnitude limit: e.g., the fraction of variance in PCi is only ~ 3% larger in the Iab < 22.5 
sample relative to the Iab < 24 sample. 

Several methods are proposed in the literature to establish the number of PCs which 
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are sufficient to fully desc ribe the proper ties of the sample; all methods require some degree 
of judgment. For example iKaiser I (jl960l ) proposes the rule-of-thumb of rejecting all compo- 
nents that contain less power than the var i ance e xpected for uncorrelated variables (in our 
case of 5 variables, less than 20%); IJoUiffd (119721 ) adopted instead a lower threshold value. 
To classify galaxy structure, we will use the ffist three PCs; these explain 92% of the total 
variance. The ZEST classification grid is therefore constructed in a three-dimensional space; 
specifically, galaxies are ranked, according to their structural properties, in unit cubes of the 
PCi, PC2, PC3 ZEST space. 



2.3.1. Morphological classification of the PC1-PC2-PC3, unit cubes 

(COSMOS) Galaxies with different structural properties occupy different regions of the 
PC1-PC2-PC3 space. For example, PCi is dominated by C, M20, and G. Highly negative 
values of PCj^ are populated by highly centrally-concentrated galaxies. PC2 is, on the 
other hand, dominated by ellipticity and asymmetry: round asymmetric objects are found 
at negative values of PC2, and symmetric flattened systems are preferentially located at high 
positive values of PC2. PC3 is also mostly a combination of asjTiimetry and ellipticity, but 
in PC3 these parameters contribute both with positive weights to the absolute value of PC3. 
Highly asymmetric and elongated objects are thus located at high values of PC3. 

To associate a (dominant) morphological class to different regions of the PC1-PC2-PC3 
space, the latter was partitioned into a regular 3-D grid with unit steps in each of the 
coordinates, and all galaxies in our COSMOS sample within each unit PC1-PC2-PCS cube 
were visually inspected. Each unit cube was then assigned a morphological type (T=l, 2 
or 3 for early-type, disk and irregular galaxies, respectively; face-on SO galaxies would of 
course be found in cubes classified as T=l, while more inclined SO galaxies would be in 
cubes classified as T=2), and a dumpiness parameter (in unit steps, ranging from = for 
smooth surface density distributions, to 2 for very clumpy morphologies). For galaxies of 
T=l or =2, we furthermore assigned an elongation parameter (in unit steps, from for face- 
on galaxies, to 3 for edge-on galaxies), and an irregularity parameter (in unit steps, ranging 
from for regular surface density distributions, to 2 for disturbed T=l or =2 morphologies). 
A measure of the galaxy sizes (i.e., their Petrosian radii, see Section rA.2p is also available 
for all ZEST-classified galaxies, as a byproduct of our analysis; see Appendix IA.2p . 



^M2o is always < 
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2.3.2. Bulge-to-disk ratios: Parametric surface brightness fits 

To refine the ZEST structural classification of T = 2 disk galaxies, we use additional 
information that we have available for a subsample of the COSMOS galaxies considered 
in this paper, i .e., th e GIM2D single-Sersic fits to the Iab < 22.5 COSMOS galaxies of 



Sargent et al. I (l2006l . this volume). In particular, we use the statistical distribution of Ser- 
sic index n within each unit cube of PCi-PC2-PC^ with a T=2 classification, to assign a 
hulgeness parameter to each T=2 cube. Specifically, the T=2 unit cubes are split in four 
bins, i.e, T= 2.3, 2.2, 2.1 and 2.0, depending on the value of the median Sersic index n of 
the galaxies in that cube (T=2.3, 2.2, 2.1, 2.0 for < Umedian < 0.75, 0.75 < Umedian < 1-25, 
1.25 < Umedian < 2.5, and rimedian > 2.5, respectively). This refinement of the Type-2 
ZEST classification grid can be interpreted to correspond to a four-bins classification of 
disk galaxies in terms of their bulge-to-disk ratios, with the T=2.0, 2.1, 2.2 and 2.3 cubes 
hosting galaxies with decreasing bulge-to-disk ratio (Type 2.0 = bulge-dominated galaxies, 
including relatively inclined SO galaxies, and Type 2.3 = bulgeless disks). 



2.4. Summary: The ZEST PC1-PC2-PC3 classification scheme and grid 

We summarize the final ZEST classification scheme in Table [21 and the COSMOS- 
calibrated ZEST grid in Table [31 

To show the performance of ZEST in disentangling galaxies with different structural 
properties, we plot in Figure [31 a representative selection of the Iab < 24 COSMOS galax- 
ies that occupy four different unit cubes of PC1-PC2-PCS. Furthermore, in Figure [H we 
show, in sequential planes of constant PC3 = —2, —1, 0, +1, +2, +3, a representative galaxy 
in each of the {PC^-)PCi-PC2 unit bins. Arrows in the bottom-left corners of each panel 
(i.e., PC3 =constant plane) indicate the directions of the steepest (positive) variation for the 
quantities specified as labels of the arrows; e.g., arrows labeled as "bulgeness" or "irregular- 
ity" respectively show the direction, across the given PC3=constant plane, of the maximum 
increase of the "degree of bulgeness" and of the irregularity of the galaxies populating that 
plane. 

In Figure[5lwe summarize the COSMOS-calibrated ZEST classification grid in a schematic 
way. In each PC3 =constant plane, different symbols represent the different morphological 
types, elongation and bulgeness parameter. As indicated in the Figure, ellipses represent the 
T = 1 early-type galaxies; concentrical double circles indicate T = 2 disk galaxies, and stars 
represent T = 3 irregular galaxies. The size of the internal ellipse of T=2 galaxies is propor- 
tional to their bulgeness parameter. The continuity of properties in the PC1-PC2-PC3 space 
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is immediately evident from Figure [51 For example, the transition from early-type morpholo- 
gies (T = 1), to bulge-dominated disks (T = 2.0), to pure-disks galaxies (T = 2.3) is smooth 
both across the PC1-PC2 planes, and along the PCs direction. The bulge-dominated galax- 
ies are found preferentially at intermediate values of PC3, and, as expected, at low values of 
PC,. 



3. The performance and reliability of ZEST 



Before demonstrating the robustness of the ZEST classification it is important to stress 
that generally, galaxy appearance depends on the rest-frame wavelength at which it is ob- 
served. Since only F814W ACS images are available for the COSMOS galaxies it is important 
that any comparison of morphology for galaxies at different redshifts is treated with care. In 
the low-to-intermediate redshift regime (0.2 < z < 1.0) that is the focus of this paper, the 
central wavelength of the F814W filter covers the rest-frame 4000-6 700A window, where mor 
phological i^— corrections are negligible for most galaxies (see, e.g. jLotz. Primack fc Madau 
20041 : ICassata et al.lboosi ). 



3.1. Testing ZEST on z = galaxies 



We assess the perfor mance and reliability of the ZEST classification grid by applying 



it to the iFrei et al.l (119961 ) sample of 8 2: = representative galax ies that (a) have Hubble 



types available from the RC3 catalog (Ide Vaucouleurs et al.lll99ll ). (b) have been observed 
(at the 1.1 meter telescope of the Lowell Observatory) in the Bj band (Ae// = 4500A) with a 
pixel scale of 1'.'35 per pixel , and (c) have PSF-FWHM smaller than S'.'O. We excluded from 
the analysis three galaxies for which the available Frei's images were too small to reliably get 
an estimate of the background. Frei's galaxies span Hubble types from Ellipticals (T = —5) 
to Sd (T = 10); they have been used 01 



Bershadv et al. 


2000; 


Simard et al. 


2002 



Lotz. Primack fc Madau Il2004f ). 



In Figure [6] we show the fraction of objects with a given RC3 classification (E, SO-Sab, 
Sb-Scd and Sd and later) that have respectively Type=l, 2.0, 2.1, 2.2, 2.3, and T=3 ZEST 
classification. The comparison between the two classifications is excellent, and highlights 
the power of the automatic ZEST classification scheme to recover the physically-motivated 
Hubble types of galaxies. 



There are only a few objects which have significantly different classifications between 
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RC3 and ZEST. In Figure [7] we show the postage stamps of these most discrepant galaxies. 
A detailed analysis of these galaxies shows that: 



disk component ( 


Micl 


lard & Marchal 


1994; 


Scorza & Bender 


1995; 


Mizuno & Oikawa 


1996; 


Emsellem et al. 


2004). 



NGC 4710 has a ZEST Type= 2.2 and a Hubble Type SO. As the image shows (see 
also, iMichard &: Marchal Ill994l ). this edge-on galaxy displays a very bright ring and 
an important equatorial dust lane, which causes its visual classification to be highly 
uncertain: indeed, in the UGC catalog this galaxy is classified as an SOa. 

The two galaxies classified as Type 1 by ZEST and with an RC3 classification of SO-Sa 
are relatively face-on galaxies. Their surface brightness is rather smooth with no visible 
spiral arms, or star formation, with the exeption of a smooth ring around NGC 4340. 

NGC 4088, classified by ZEST as Type= 3 and by the RC3 as Sb-Scd, is described as 
an irregular and distorted spiral by Dahari (1985). 

Finally, for the few galaxies which are classified by ZEST as bulge dominated (Type 2.0) 
galaxies and have RC3 classification of Sb-Scd, we performed a single Sersic fit to their 
surface brightness profile. They have concentrated light distributions, with Sersic's 
index n > 3, which confirms our ZEST classification as Type 2.0 galaxies. 



3.2. Error bars in the ZEST morphological classification of COSMOS galaxies 

In Appendix |B] we discuss in detail the uncertainties and the systematic errors in the 
measured structural parameters as a function of signal-to-noise ratio (S/N). To do so we use 
a sample of bright COSMOS galaxies, which are progressively dimmed to fainter magnitudes 
(lower S/N). Here we summarize to what extent the COSMOS-calibrated ZEST morpholog- 
ical classification grid is affected by the S/N of the individual galaxies. We use the same 
sample of bright and progressively S/N-degraded test galaxies described in Appendix |Bl and 
compute for each of them all the parameters involved in the ZEST classification, both on 
the original and on the progressively S/N-degraded images. The PCi {i = 1, 2, 3) values are 
then computed for each galaxy at each S/N level, and the ZEST morphological classification 
corresponding to the relevant unit cube of PC1-PC2-PC3 is assigned to each of the original 
and artificially fainted galaxies. The change in galaxy type (AT) that occurs due to degraded 
S/N is finally computed as a function of magnitude. 
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The ZEST classification is robust down to Iab = 24. For magnitudes Iab < 22.5, 
more than 90% of galaxies do not change morphological class. The remaining few percent of 
galaxies changes morphological type by smoothly moving through the PC1-PC2-PC3 space; 
indeed, the change of ZEST morphological Type with varying S/N happens typically for 
galaxies which are originally classified in PC1-PC2-PC3 cubes that are close to a "type 
transition-wall" in the PC's space. The fraction of galaxies with AT = remains larger 
than ~ 75% down to magnitudes Iab = 23.0, and even in the highest magnitude bin (23.5 < 
/ < 24), the fraction of galaxies with AT = remains of order 65%. This is illustrated in 
Figure [U where we show the distribution of the average absolute variation in PCi, namely 
(A(PC)) = (^. \PCij - PCi,o|)/3 (with i = 1,2,3, PQ^q the initial values of the PC^'s, 
and PCi J the measured PCis after S/N degradation). Split in four magnitude bins (from 
Iab = 22, top-left panel, down to Iab = 24, bottom-right panel), the solid histograms show 
the (A(PC)) distribution for all galaxies in the considered magnitude bin, and the hatched 
histograms show the (A (PC)) distribution for galaxies which change morphological class 
due to S/N degradation. 

Figure [9] shows a few examples of how galaxies move through PC space as the S/N of the 
galaxy images decreases. Initial and final values for the PCi coordinates of the example test 
galaxies are indicated at the beginning and end of the track that describes the movement, in 
the PC space of the specific test-galaxy. Note that galaxies can have a (A (PC)) as high as 
(A (PC)) = 2, without changing their morphological class (as their path in PCi space occurs 
within a region uniformly classified with a specific morphological type). As stressed above, 
virtually all galaxies that change morphological class are located, to start with, in bins of 
PC1-PC2-PC3 that border with bins with a different morphological classification. 

The most noticeable effect of S/N degradation is a contamination, at the faintest mag- 
nitudes, of at most ~ 30% from low S/N ear l y type galaxies to the Type 2 o r Type 3 galaxy 



populations (see also, [Abraham et al.l 119961 : iLotz. Primack fc Madau 1 120041 ). This suggests 



that at most 30% of early-type galaxies at the faintest magnitudes could drop from the 
early-type sample and be misclassified as disk or irregular galaxies due to their lower S/N 
values. We discuss in Section 15.51 the implications of this effect on the redshift evolution up 
to z ~ 1 of the LFs of COSMOS galaxies structurally classified with ZEST. 



4. The advantage of the ZEST classification scheme over other approaches 



Popular classifications of (nearby and high— z) galaxies in the past few years have been 
typically based on co nibinations of two or three of th e non-parametric diag nostics that are 
used in ZEST (e.g., iLotz. Primack &: Madau I l2004l : iFerreras et al. I l2005l . and references 
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therein) or on a threshold in Sersic index n (most noticeably the SDSS defined early-type 
galaxies the objects with C > 2.87 oi n> 2.5, and late-type galaxies all other galaxies). Such 
approaches can be seen as "lower-dimensionality" projections of the ZEST grid, and lead 
to galaxy samples that are affected by large contamination of systems with rather different 
structural properties. 

As an example, in Figure [10] we show the two-dimensional planes defined by G and 
M20 (top row), C and M20 (middle row), and C and A (bottom row). All plots show the 
total density of galaxies in grey scale. In the different panel we indicate the region of space 
that is populated by a given ZEST morphological class; in particular, the color contours 
enclose ~ 99% of the COSMOS galaxies with the specified ZEST morphological class. To 
avoid crowding we show in the left panels the location of the T = 1 early-type galaxies (red 
contours), in the central panels the T = 2 disk galaxies (blue contours, with decreasing shade 
of blue from the bulge-dominated T = 2.0 galaxies to the T = 2.3 bulgeless disks), and in 
the right panels the T = 3 irregular galaxies (green contours). 

It is clear from Figure [TD] that there is a high level of contamination by different galaxy 
types in all regions of these two-dimensional planes. For example, on the C — M20 plane, 
T = 1 early-type galaxies form a tight sequence, well visible in the left panel of Figure [T0| 
however, such a sequence includes (not only the majority of T = 2.0 bulge-dominated disk 
galaxies, but also) a substantial fraction of lower bulge-to-disk ratio (i.e., T = 2.2,2.1) 
galaxies. For example, if the galaxies that lie above the solid black line on the C — M20 
plane were classified as "early-types", ZEST would return the following galaxy population 
mixture in the selected region: 38% of desired T = 1 early-type galaxies, 23% of T = 2.0 
galaxies, and 26% and 10% of T = 2.1 and of T = 2.2 galaxies, respectively. Even summing 
up together the T=l and 2.0 galaxies (which might be desired for some science applications), 
still the selected sample would be affected by a contamination of order 40% contributed by 
lower "bulge-to-disk ratio" galaxies. Similar levels of contamination are found when the 
other two-dimensional planes of Figure [10] or a simple cut in Sersic index n are used to 
morphologically-classify galaxies, as it has been done in most of the previous literature. 

5. A first application of ZEST: The evolution since 2; ~ 1 of the LFs of 
morphological early-type, disk and irregular COSMOS galaxies 

We derive the rest-frame 5— band LFs of ZEST-classified COSMOS galaxies brighter 
than / = 24 (see Appendix lA.ll for details on the COSMOS sample). Given the large 
number statistics of (the fraction of survey area that we are considering for) COSMOS, we 
can compute the LF $(M, z, T) of each morphological class in four different redshift intervals. 
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We use the ZEBRA Maximum Likelihood photometric redshifts of Feldmann et al. 
(2006) to derive the COSMOS LFs. These photometric redshifts have an accuracy of Az/(1 + 
z) ~ 0.03 in comp arison with the zCOSMOS spectroscopic redshifts of Iab ^ 22.5 galaxies 



f lLillv et al.ll2006f ): the accuracy of the ZEBRA redshifts degrades to Az/{1 + z) ~ 0.06 
down to our magnitude limit of Iab = 24 (see Appendix IA.3I for details). The application 
of ZEBRA to our sample shows some dependence of the resulting photometric redshifts on 
whether small (of order ~ 0.05 magnitudes or smaller) systematic offsets, which are detected 
by the code, are applied to the photometric calibration of the COSMOS Subaru data. This 
has however no substantial impact on our resulting LFs, which we present below as computed 
both with and without corrections for these photometric offsets. 

The final Iab < 24 COSMOS sample that we study below consists of 30760 galaxies 
classified with ZEST as T = 1 (2497 objects), T = 2 (26873 objects), and T = 3 (1390 
galaxies), which have photometric redshifts in the range 0.2 < z < 1.0. 



5.1. Luminosity Functions: Definitions 

We es timated the rest-frame--B ga laxy LF in different redshift bins, using the 1/Kn 



estimator (jSchmidt Ill968l : iFelten Ill976l ). 



According to the original 1/Knax formalism, the number of galaxies per unit comoving 
volume in the range of absolute magnitudes dM, at redshift z, and morphological class T 
can be written as: 



j $(M,2,r)dM = J2 



V 



where the sum is over all galaxies in the specific range of redshift, absolute magnitude, and 
morphological class. Knax,i is the maximum comoving volume within which the galaxy i 
could still be detected according to the apparent magnitude limits of the survey, which, in 
our case, is given by 16 < / < 24. The Kiax,i is computed for each galaxy according to: 

l-mm{zu ,Z24) 

K^ax,i = ^ / -r-dz, (9) 

where zu and zl are the upper and lower redshift of the considered redshift bin, and 2:24 and 
ziQ are the redshifts at which a galaxy of a given rest-frame 5— band absolute magnitude 
and a given SED would have /—band apparent magnitude of 24 and 16, respectively. The 
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values of 2:24 and ziq were computed, for each galaxy, taking into account the A;— correction 
resulting from the best fit SED provided by the ZEBRA fit. No evolutionary correction was 
applied. Q is the effective area of the survey, corresponding to 0.74 square degrees. The 



Since we are working with a magnitude-selected rather than a volume-limited sample, 
the highest and lowest 5— band magnitudes at each redshift depend on the galaxy SED. 
This is evident in Figure [19] of Appendix IA.41 where we show the color K-correction (i.e., 
difference between the observed /— and the 5— band rest-frame magnitude), as a function 
of redshift, for templates of different photometric type (color-coded from red for early type 
galaxies up to blue for starbursting galaxies). The effect vanishes at redshift z = 0.8, where 
the central wavelength of the F814W filter matches exactly the rest-frame i?-band. Since 
correcting for this bias needs an a priori assumption on the color distribution of the non- 
detected galaxies, we limit the computation of the LP in each redshift bin to the luminosity 
range for which we are complete, regardless of galaxy colors. The used magnitude ranges 
are listed in Table HI 



About 2.4% of the ACS-detected galaxies do not have a match in the ground-based 
catalog; therefore, no information on their spectral energy distributions - and thus photo- 
metric redshifts - can be derived. Furthermore, about 2% of the ACS-detected sources are in 
close pairs, and, particularly, are at a distance from each other smaller than 0'.'6; therefore, 
these sources are detected as a single object in the ground-based images (see Appendix lA.ll 
for details). The photometric redshifts estimated for these objects are thus unreliable, even 
for those pairs which are really physically associated and thus at the same distance, since 
the ground-based images mix the light from both galaxies, which could have very different 
stellar population. Therefore, these galaxies in close pairs also need to be excluded from the 
analysis of the COSMOS LFs. 

In Figure [TT] we show the /—band magnitude distribution of the ACS-galaxies with no 
detection in the ground-based catalog (top panel), and the magnitude distribution of the 
ACS-galaxies for which the same groundbased identification was associated with multiple 
ACS sources (bottom panel). Both distributions are normalized to the number of galaxies 
in the total sample in each magnitude bin. The fraction of missed objects stays relatively 
constant down to a /—band of ~ 22, but then increases toward fainter magnitudes. Although 
the fraction of missed objects is small, its dependence on the observed /—band magnitude 
may introduce biases in the computation of the LF. 



Poissonian errors from galaxy counts on $(M, z, T) are given by 0"$ 




5.2. Corrections for missing sources 
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To assess the potential impact of such biases, we therefore computed two versions of the 
LFs, respectively assuming that: 

- All galaxies with no redshift are outside the 0.2 < z < 1.0 range. We refer to these 
LFs as to the uncorrected LFs. 

- The missing galaxies have the same redshift distribution than the galaxies with avail- 
able redshift. In this case we can calculate the statistical weight (xi) needed to correct 
the Vmax value of galaxies with known redshift. We calculat e y,; as a function of th e ob- 



served J— band magnitude, following the approach used by lWillmer et al. I ( l2006l ). We 
consider magnitude bins of 0.5 magnitudes, and for each magnitude bin we compute: 
(z) The number of galaxies with photometric redshift in the considered redshift range 
{Nz, with 0.2 < z < 1.0); (ii) The number of galaxies with z > Zh {Nzh, with z^ = 1.0); 
and [iii) The number of galaxies with z > zi [Nzi, with zi = 0.2). Under the assump- 
tion that objects without a redshift estimate have the same redshift distribution of 
the entire sample, the probability for a galaxy of magnitude of being in the redshift 
range 0.2 < z < 1.0 is given by the number of galaxies with good redshift estimates 
in that range, divided by the sum of the number of objects with good redshift plus 
the number of objects with redshift both lower and higher than the considered limits, 
i.e.. Pi = Nz/{Nz + Nzh + Nzi). Finally, the weight Xi of each galaxy with estimated 
redshift is the sum of the Pi of all galaxies in the relevant magnitude bin, divided by 
the number of galaxies in that bin with 0.2 < z < 1.0 (i.e., ^Pi/N^). The Knax,i 
that is associated with each galaxy, weighted to account for objects with no available 
redshift, is simply Xi/'^max/ii and the Poissonian errors on the LF, in each magnitude 
bin, are given by af = XlXi/Knaxi- following we will refer to the so-obtained 

LFs as to the corrected LFs. 



In Figure [12] we show the LFs of the total sample of Iab < 24 COSMOS galaxies, as 
well as the LFs for the different ZEST types. Specifically, we show the comparison between 
the uncorrected (dotted line) and the corrected (stars) LFs. The figure shows, from top to 
bottom, the LF integrated over all morphological types, and the LFs for the T =1 (early- 
type), T =2 (disk) and T =3 (irregular) galaxies, respectively. Each column shows the LFs 
calculated in different redshift bin: 0.2 < z < 0.4 (first column), 0.4 < z < 0.6 (second 
column), 0.6 < z < 0.8 (third column), and 0.8 < z < 1.0 (last column). The plotted error 
bars show the Poisson errors only. Uncertainties due to sampling variance are not included 
in the figures; these are discussed in Section 15. 4[ In each redshift bin, the LFs are shown 
down to the absolute 5— band magnitude at which the sample is complete, regardless of 
galaxy colors, so that both red and blue galaxies are sampled in an unbiased way at each 
redshift, down to the faintest magnitude bin. 
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5.3. Impact on the LFs of the uncertainties in the photometric redshifts 

In order to check to which extent the limited redshift accuracy affects the derived LFs, we 



perfor med a set of simulations using the COSMOS mock galaxy catalogs of lKitzbichler et al. 



(120061 ): these are designed to reproduce in detail the COSMOS survey. From these mock 
catalogs, we first extracted all galaxies with observed magnitude Iab < 24, covering an 
area on the sky equal to that analyzed in this paper, and then generated 100 simulations of 
the galaxy catalog, perturbing each redshift with an error derived by randomly sampling a 
Gaussian with a a equal to the one of the photometric redshift uncertainty in the relevant 
redshift bin. We considered ^"^/(l + 2;) ~ 0.03 for galaxies brighter than Iab = 22.5, and 
az/{l + z) ~ 0.06 for fainter objects (see Appendix lA. 30 . The LFs of the original catalog and 
of the 100 re-simulated samples were then calculated using the same procedure described in 
Section 15.11 The results of these tests are presented in Figure [131 where, in each panel, the 
solid line with Poissonian error bars represents the LF computed using the "true" redshifts, 
and the open diamonds represent the median of the 100 realizations. The shaded grey area 
associated with the simulated volume densities represent the 16*^* and 84*^^ percentiles of 
the simulated distributions within each magnitude bin. Figure [13] clearly shows that the 
dominant effect of the redshift uncertainty are the systematic trend of (a) populating the 
bright end of the LF, where few or no galaxies are present in the original sample, and (6) 
slightly underestimating the density of galaxies around the knee of the LF (i.e., around M^,). 
However, even at the bright-end of the LF, where the effect is the strongest, the real and 
the simulated LFs are well within 2cr$, implying that the uncertainty in the photometric 
redshifts does not affect our conclusions on the evolution of the LFs. The largest source of 
uncertainty, especially in the lowest redshift bin, is due the small volume sampled by the bin, 
where the contribution from large scale structure variation are significant. This explains the 
relative large difference in the LFs calculated from two independent mock catalogs (shown 
with dashed lines in Figure [T3!l . 

In Figure [12] we show, as stars, the corrected LFs derived adopting the ZEBRA photo-z 
calculated after correcting the photometric catalogue for the small systematic offset, and, 
shaded in grey, the uncertainty on the LFs that arises when using the ZEBRA photo-z esti- 
mates obtained without correcting the photometric catalogue (see Appendix IA.3p . Although 
the LFs estimated with the two sets of photo-z's differ somewhat (see Table [1] for quanti- 
tative differences in the corresponding Schechter fits), the differences do not substantially 
impact our conclusions. Furthermore, Figure [12] illustrates that the LFs derived with the 
two approaches described in Section [5^ are in very good agreement, demonstrating that the 
results are independent of the (small) fraction of objects excluded for not having a photo- 
metric redshift estimate. In the following we therefore focus our discussion entirely on the 
corrected LFs (stars in Figure [T^j). 
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5.4. Results: The evolution of the LFs for the different morphological classes 



The first column of Figure [12] shows that, in the lowest redshift bin, the shape of the LP 
varies significantly for different morphologically-selected galaxy samples. In particular, the 
LF of the global galaxy sample and the LP of T = 2 disk galaxies keep increasing toward 
faint magnitudes, and have a rather similar behavior. On the other hand, the LFs of the 
T = 1 early-type galaxy population, and of the T = 3 irregular galaxies appear, within the 
errors, to almost flatten for magnitudes Mb > —19.5. 

We note that our COSMOS LFs computed in the flrst redshift bin {z = 0.2 — 0.4) are 
rather susceptible to the effect of large scale structures, as the size of the studied COSMOS 
fleld at z = 0.3 is only ~20 comoving Mpc. To quantify this effect we compare the LF 
for the global galaxy population with the LFs derived in the same redshift bin by other 
surveys. A flt with a Schechter function to the COSMOS LF in this redshift bin gives 
the following best flt parameters: = 0.005 ± 0.001Mpc-^ Mb,.^ = -21.0 ± 0.2 and 
a = —1.26 ± 0.15. The fa i nt-en d slope a is in good agreement with the value recently 

-1.3) on the basis of the DEEP2 spectroscopic 



publish ed bviwmmer et ai. I ll^uubi . awiUmer 

survey (jPavis et al.ll2003l ). The DEEP2 LF is based on spectroscopic redshifts; therefore, the 
consistency of the DEEP2 LFs with the one that we have presented highlights the reliability 
of our photometric redshift estimates. In the 0.2 < z < 0.4 bin, the comparison between 
our value of and those derived from other surveys demonstrates that the average density 
in the C OSMOS fleld is systena atically higher than in other studies at sir ailar redshif t s. For 
example, Iwillmer et aD J2006h flnds 0* = 26.39^^;^^ x 10"^ Mpc-^ and IWolf et aP J2003 ) 
flnds (j)^ = 18.27 ± 12.98 x 10"^ Mpc'^ in this redshift regime, for the C0MB017 survej- 
The overdensity in the COSMOS fleld at 0.2 < z < 0.4 is conflrmed by the spectroscopic 
(Lilly et al. 2006) and photometric (Feldmann et al. 2006) redshift distributions. Similarly, 
a large overdensity aX z = 0.75 is present in the COSMOS fleld (Scoville et al. 2006), which 
is also detected by XMM ( iHasinger et al. II2006I . this volume). 



Figure [T^ shows that the LF of the global galaxy population does not change signiflcantly 
as a function of redshift. Qualitatively, the behavior of the global LF can be described by a 
brightening of the global galaxy population from redshift z = 0.3 to z = 0.9, with the volume 
density of the entire galaxy population staying basically constant from redshift z = 0.5 to 
z = 0.9. The Type 2 disk galaxy population follows closely the behavior of the global sample 
as a function of redshift, and can be described with a rather constant number density, and 
a brightening of the stellar populations. 



^The C0MB017 survey covers a total area of 0.78 deg^, in t hree fields. The m easured <&* in the different 
fields ranges from 13.03 x IQ-* Mpc-^ to 21.17 x lO"'^ Mpc-^ |Wolf et al.ll2003l ) 
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In Figure [m we show the corrected LFs of T = 2 disk galaxies sphtted by their bulgeness 
parameter. As in Figure [121 different columns show the results for different redshift bins from 
0.2 < z < 0.4 to 0.8 < z < 1. The LFs of different types of disk galaxies are plotted in 
different rows, from the bulgeless T = 2.3 disks in the top row, to the bulge-dominated 
T = 2.0 disk galaxies in the bottom row. In the lowest redshift bin, as the contribution 
of the bulge component increases (i.e., moving from top to bottom in Figure [T^ . the faint 
end of the LF becomes flatter, and the LF becomes more similar to that of the Type 1 
galaxies. At low redshifts, the bulge-dominated disk galaxies contribute most of the light 
at magnitudes brighter than Mb = —21.5, while their contribution is strongly reduced at 
magnitudes fainter than Mb = —19.0. The LF of T = 2.0 bulge-dominated galaxies is very 
similar in both shape and normalization to the LF of the T = 1 early-type galaxies. 

5.5. Anchoring the COSMOS evolution to the local universe: Comparison 

with a complete SDSS-based sample 

The presence of the z < 0.4 overdensity in the COSMOS field hampers the interpretation 
of the evolution of the LFs over the studied redshift range; the small volume sampled by the 
COSMOS survey at redshifts z < 0.2 is furthermore inadequate to constrain the LFs at such 
late epochs. To correctly compare the COSMOS data with the z = Universe, we used the 
SDSS-based volume-limited sample of 2; = galaxies (appropriately "redshifted" to z = 0.7, 
hereafter indicated as SDSS2=o.7) constructed and extensively discussed by Kampczyk et al. 
(2006); we briefly summarize its properties in Appendix [Cl 

We analyzed the SDSS2=o.7-simulated images following the identical procedure that 
we adopted for the COSMOS galaxies: We classified the SDSS galaxies according to their 
position in the 3-dimensional PCi space, and then computed the LFs of the morphologically- 
selected SDSS2=o.7 sample, as described in Section [STTl This SDSS2=o.7 sample provides the 
properties that local galaxies would show, at redshift z = 0.7, assuming no evolution between 
z = 0.7 and z = 0. We therefore stress that applying the identical classification and LF- 
derivation procedures to both our COSMOS and our SDSS-based samples guarantees that 
any systematics present in the former is also present in the latter; thus, the comparison 
between the two samples allows us to discover evolutionary trends that are unaffected by 
such systematics. 

The comparison of the COSMOS 0.6 < z < 0.8 LFs (stars) with the SDSS2=o.7 samples 
(dash-dotted curves) is shown in Figure [151 The comparison shows that pure-luminosity 
evolution could explain the differences between the global local SDSS galaxy population 
and the global COSMOS population at 2; = 0.7. To quantify this statement, we compute 
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the brightening of the global local galaxy population that is required to overlap the global 
SDSS^=o.7 LF with the global COSMOS LF at redshift z = 0.7, which is of order 0.95 ± 0.1 
magnitudes. This amount of brightening is applied to the SDSS2=o.7 LFs in the various plots 
of Figure \TE\ (dashed curves). The agreement between the brightened SDSSz=o.7 LF derived 
for the global galaxy population and the global COSMOS z = 0.7 LF is remarkable. Thus, 
the evolution of the global galaxy population since z = 0.7 can be interpreted as due to only 
a brightening of the stellar populations, with no evolution in the average volume density of 



(2006 


) and also 


Faber et al. 


(2005) 



of the DEEP2 spectroscopic survey. 

The Schechter fits to the total. Type 1, Type 2 and Type 3 COSMOS LFs in the 
0.6 < z < 0.8 redshift bins are reported in Table [Hand, as solid black lines, also in Figure [T^ 
In particular, the Schechter fits to the z = 0.7 bin are also plotted in the lower redshift bins, 
after applying the amount of luminosity evolution that is required to overlap the 0.6 < 
z < 0.8 COSMOS LFs with the local SDSS LF (which corresponds to a 5— band evolution 
proportional to ~ 1.3 x 2;, as in Figure [T5l) . Despite the noticeable effect in the lowest 
redshift bin of the known structure in the central part of the COSMOS field, the shapes of 
the COSMOS LFs at z ~ 0.7 well describe, for all the different morphological populations, 
the LFs at lower and higher redshifts. 

The further inspection of Figure [15] shows however, qualitatively, some interesting evo- 
lutionary trends in the individual morphologically-selected galaxy samples. In particular: 



- Not surprisingly (since they dominate the global galaxy population), the evolution of 
the LF of the Type 2 disk galaxy population is very similar to that observed in for 
the global population: the LF of disk galaxies is consistent with a pure luminosity 
evolution of 0.95 magnitudes up to z = 0.7, with a constant average volume density. 

- Within the errors, the bright-end {Mb < —21.5) of the LF of morphologically-classified 
Type 3 irregular galaxies aX z = 0.7 is also consistent with that of the local SDSS 
irregular galaxies (brightened by a similar amount as above). At fainter luminosities, 
Mb > —20.5 we detect however an excess of a factor of order ~ 3 in the number density 
of faint irregular galaxies at z ~ 0.7 relative to the local universe. 

- At bright magnitudes Mb < —21.5, there is also a good agreement between the LF of 
z = 0.7 COSMOS Type 1 early-type galaxies and the LF of similarly-classified SDSS 
galaxies brightened by 0.95 magnitudes. This evolution is similar to that expected the- 
oretically from the passive evolution of a coeval stellar population formed well before 
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2; ~ 0.7, and also to the evolution that is observed from studies of the surface bright- 
nesses of early-type galaxies (Rigler & Lilly 1994, Schade et al. 1999, Mcintosh et al. 
2005, Holden et al. 2005). This result supports a scenario where most of the massive 
early-type galaxies are in place at z = 0.7, with the caveat that we cannot be sure 
that all of the z ~ 0.7 population are actually the progenitors of early-type galaxies 
(e.g., not all of the T=l galaxies may be consistent with evolving onto the Kormendy 
relation aX z = 0,Ferreras et al. 2005; we investigate this issue further in Scarlata et al. 
2006b). At fainter magnitudes, Mb > —20.5, however, there is a drop of a factor ~ 2 
between z = and z = 0.7 in the number density of morphologically-classified early- 
type galaxies. This result is potentially in contrast with a significant contribution of 
elliptical-elliptical mergers to build the most massive early-type galaxies at relatively 
recent epochs, and favors a down-sizing scenario for early-type galaxy formation, in 
which the most massive ellipticals form and dynamically-relax at earlier epochs than 
the lower-mass ellipticals. This is also suggested by some studies of ellipticals in the 
local universe (e.g., CaroUo et al. 1993; Thomas & Davies 2006) and at high redshifts 
(e.g., Treu et al. 2005, Daddi et al. 2005). 

The ZEST misclassification at the faintest magnitudes discussed in Section 13.21 is not 
responsible for the observed deficit of Type 1 early-type galaxies and increase in the number 
density of Type 3 galaxies, since: (1.) at most ~ 30% of Type 1 galaxies could be misclassified 
as Type 2 or Type 3 galaxies; this fraction is too small to explain (a) the missing factor of 
two of Type 1 galaxies at redshift z = 0.7; and (b) the factor of 3 excess in Type 3 galaxies 
observed in the same redshift bin. (2.) The SDSSj;=o.7 galaxy sample accurately simulates 
the COSMOS ACS images. Therefore, the same fraction of COSMOS and SDSS2=o.7 galaxies 
would be misclassified by ZEST due to low S/N. Thus, the relative comparison of the two 
samples is not affected by misclassification due to low S/N. 

We postpone to future publications more sophisticated analyses of the LFs of COSMOS 
galaxy samples that are based on the structural/morphological classification obtained with 
ZEST. 



6. Summary 

The Zurich Estimator of Structural Types (ZEST) uses the Sersic index n of the fits 
to the galaxy surface brightness distributions, and five basic non-parametric diagnostics to 
quantify the properties of galaxy structure. The novelty of ZEST is to use a Principal 
Component Analysis to retain the full information provided by the entire set of diagnostics 
while reducing to three the dimensionality of the structural parameter space (of axes PCi, 
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PC2, and PCs). The ZEST scheme morphologically classifies galaxies in three main types 
(early-types, disk galaxies and irregular galaxies), assigns a bulgeness parameter related to 
the B/D ratio (quantized in four bins) to disk galaxies, and ranks galaxies according to the 
elongation, irregularity and clumpyncss of their light distribution. Most important to stress is 
that, in contrast with other approaches to morphological classifications that utilize only two 
or three of the basic structural parameters that enter the ZEST PCA, the ZEST PC1-PC2- 
PC2 classification scheme substantially reduces the contamination to each morphological 
class from different galaxy types. We have used a sample of ~ 56000 Iab < 24 galaxies in 
the central 0.74 deg^ of the COSMOS field to cahbrate the ZEST classification grid, and, so 
far, ascertained the robustness of the ZEST classification down to the faintest magnitudes 
of our sample; work to stabilize the results down to even fainter magnitudes is in progress. 

As a first application, we use our ZEST-classified sample of COSMOS galaxies to study 
the evolution since 2; = 1 of the Luminosity Function of the early-type, disk (of different 
bulgeness parameter) and irregular galaxies. Our analysis shows that the average volume 
density of disk galaxies remains constant through the studied redshift range, although the 
stellar populations of these systems are brightened at earlier epochs. In contrast, only 
the bright, Mb < —21.5 end of both the irregular and the early-type galaxies remains 
roughly consistent with the brightened LP of local galaxies. At faint magnitudes, both 
morphologically-classified irregular and early-type galaxies show substantial evolutionary 
effects from z — to z — 0.7 and above. In particular, there is an excess of a factor of order 
~ 3 of faint morphologically-classified irregular galaxies, and a deficit of a factor of about 2 
of faint {Mb > —20.5) morphologically-classified early-type galaxies at 2; ~ 0.7 relative to 
the local universe. 

The ZEST classification of our sample of 56000 COSMOS galaxies has already been 
used to study the number density evolution of disks of different sizes (Sargent et al. 2006), 
the LF-evolution of plausible progenitors of eUiptical galaxies (Scarlata et al. 2006b), and the 
evolution of the merger rate out to ^ ~ 1 (Kampczyk et al. 2006). The ZEST approach will 
be used to automatically classify the whole 2deg^ of COSMOS down to Iab — 24 (CaroUo 
et al. 2006). 
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A. Details on the COSMOS galaxies 
A.l. The sample 



A detailed discussion of the ACS data processing is given in iKoekemoer et al. I (120061 ): 



the description of t he ge neration of the ACS catalog that we use in this paper is given in 



f version 2.4.3 



Leauthaud et al. I (I20061). Briefly, the ACS catalog was produced by applying SExtractor 



Bertin fc Arnoutd 119961 ) to the reduced ACS images. Edge regions in each 
image were removed with no loss of data, since the overlap between neighboring ACS tiles 
was larger than the excluded edge regions. Two SExtractor runs were performed in order to 
avoid both blending of multiple independent sources and de-blending of large galaxies, i.e., 
a "cold" run with a configuration optimized for the detection of large, bright objects, and a 
"hot" run with a configuration optimized for faint, small objects. The two resulting samples 
were then merged together to produce a final catalog by retaining all the "cold" detections 
plus the "hot" detections which fell outside the SExtractor segmentation map of any galaxy 
detected in the "cold" run. 

We visually inspected the catalog to remove any residual over-deblending of large galax- 
ies and also false detections. Specifically, we cross-correlated the ACS catalog with the 
COSMOS CFHT /-band catalog down to / = 24.5 using the celestial coordinates of the 
objects (5(7 magnitude limit of /cfht = 24.9; Capack et al. (2006))@, and fiagged for visual 
inspection the ACS-detected objects that: (i) Had no match in the available ground-based 
COSMOS catalog; (ii) Were detected in the ground-based catalog, but had a > O.bmag 
difference between the ground-based and HST I magnitude estimates; and {in) Were linked 
with a CFHT source that had already been associated with another ACS source. About 4.5% 
of the entries in the original ACS-based catalog were removed after the visual inspection. 

Stars were then removed from the ACS-based catalog using the SExtractor CLASS_STAR 
parameter, and particularly removing sources with CLASS_STAR> 0.6. In Figure [16] we 
show the SExtractor estimate of the radius Rtot including 100% of the light of an object, 
as a function of the SExtractor ACS /—band magnitude. Objects with CLASS_STAR> 0.6 
are identified as filled circles. These form a tight, well identifiable sequence at all magni- 
tudes. The final, cleaned, ACS-based catalog that we classify in this paper includes ~ 55651 
galaxies down to Iab = 24. 



■^For each ACS source we searched for CFHT-detected galaxies within a circle of 0'.'6 radius. If more than 
one CFHT source was found within this radius, the closest CFHT galaxy was assigned to the ACS source 
as its ground-based counterpart. 
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A. 2. Data analysis 

The steps of data analysis that we performed on the 56000 COSMOS galaxies of our 
sample are summarized below. 

- For each ACS-detected galaxy, we removed in the ACS images the contamination from 
nearby objects that could affect the computation of the non-parametric coefficients. 
To create a "cleaned" version of each galaxy image, a stamp centered around the ACS- 
detected galaxy was extracted from the original ACS tile, with size equal to 3 x i?tot- 
SExtractor was used to identify the individual sources in the stamp, and the pixels 
belonging to contaminant galaxies were replaced by background noise that matched 
the properties of the noise in the original ACS image. 



The Petrosian radius (jPetrosian I Il976l . Rp) of each galaxy was calculated from the 



surface brightness profile measured from its "cleaned" stamp0. The surface brightness 
profiles were derived using elliptical annuli centered on the galaxies, in order to mini- 
mize the contribution by background pixels. The position angle and ellipticity of the 
annuli were kept fixed a t all radii, to the values measured from the image moments 



(iBertin fc ArnoutsI Il996l ). The annuli were uniformly spaced in log(-R) to increase 
the S/N in the external region of galaxies. Each galaxy profile was sampled by at 
least 10 points. We calculated Rp for all galaxies with an Half- Width-Half-Maximum 
(HWHM, as estimated by SEXtractor) larger than 1.5 pixels, so as to ensure that 
enough pixels were available for the measurement of the parameters. This cut excludes 
~ 0.2% of the entire sample of COSMOS galaxies under c onsideration. In Figur e [T71 



we show, for the sample of Iab < 22.5 galaxies analyzed by lSargent et al. I (l2006l . this 



volume), the comparison between the Petrosian radius derived in our an alysis and the 



half-light radius derived by Sargent and collaborators from GIM2D fits (ISimardlll998l ) 
to the galaxy surface brightness distributions. The solid and dashed lines represent 
the relation between Rp and R1/2 expected for a Sersic profile with n = 1 and n = 2, 
respectively. The very good correlation between the two measured quantities indicates 
that our Petrosian radii reliably reproduce the galaxy sizes, and offer the advantage 
of expanding the COSMOS galaxy sample with computed robust size measurements 
down to Iab = 24. 



"^Rp is defined as the radius at which the ratio 77 = fi{R)/{fi{< R)) between the surface brightness at 
that radius and the surface brightness averaged within that radius is equal to a given number. For a surface 
brightness distribution described by a de Vaucouleurs or an exponential profile, a value of ry = 0.2 is reached 
at Rp ^ 1.8i?i /9 and Rp ^ 2.2i?i/2j respectively (with i?i/2 the half-light or "effective" radius of the galaxy, 
Graham fc Driver 1 12005[ )). 
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• The pixels that are associated with a galaxy were defined as all the pixels within 
elliptical apertures of semi-major axis equal to the Petrosian radius of the galaxy. 
With this definition, the same fraction of light is considered, for any given shape of 
the surface brightness profile, at different redshifts, therefore allowing to consistently 
compare measurements of galaxy populations at different epochs. 

In Figure [18] we illustrate with a few examples the procedure outlined so far. In particu- 
lar, we show four COSMOS galaxies with increasing magnitudes, from I^b = 18.8 (top row) 
to Iab = 22.3 (bottom row). For each galaxy we show the original stamp (first column), 
the cleaned stamp (second column), the segmentation map (third column), and, in the last 
column, the surface brightness profile (top panel) and the ratio rj = fi{R)/ (yu(< R)) (bottom 
panel) as a function of radius R in arcseconds (solid circles). The solid horizontal line in the 
1] — R plots represents rj = 0.2. For the galaxy in the bottom row, rj decreases to a minimum 
value that is however larger than rj = 0.2, and then increses again at larger radii. This is due 
to the contamination from a nearby source, still visible even in the cleaned stamp, which 
adds fiux in the outermost isophotes (as easily seen in the /i radial profile). To obtain reliable 
Petrosian radii, we corrected these cases by fitting a linear relation to the fi — R profile in 
the outer radial bins (dashed line), and then subtracting a constant value to the surface 
brightness profile. The open circles in the r] — R plot of the bottom-row galaxy in Figure [18] 
show the Tj profile obtained after the application of this procedure, and demonstrate that we 
recover robust values of Rp even in these potentially troublesome cases. 



A. 3. Photometric redshifts for COSMOS galaxies 

Ancillary ground-based photometry in several passbands is available for the COSMOS 
field from different facilities (see summary in Scoville et al. 2006, this volume). This allows the 
estimate of accurate photometric redshifts for the COSMOS galaxies. Due to the differences 
in adopted catalogues, and the intrinsic uncertainties in photometric redshift estimates, the 
COSMOS collaboration has decided to adopt several approaches to determine photometric 
redshifts for COSMOS galaxies. Mobasher et al. (2006, this volume) present a comparison 
among several estimates for the COSMOS photo— 2;, obtained by different groups within the 
COSMOS team. 

In our work, we adopt the photometric redshifts for the COSMOS galaxies of our ACS- 
selected Iab < 24 galaxy sample by Feldmann et al. (2006). These are derived with the 
Zurich Extragalactic Bayesian Redshift Analyzer (ZEBRA). ZEBRA produces two sepa- 
rate estimates for the photometric redshifts of individual galaxies: A Maximum Likelihood 
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estimate and a 2-D Bayesian estimate. In this paper we use the Maximum Likehhood photo- 
metric redshifts. These were derived using: • Subaru photometry in six broad band filters B, 
V, G, r', i', and z' {5a magnitude hmit of ~27 for point sources in all bands; see Taniguchi 
et al. 2006, this volume, for details); • u* photometry from the MegaCam at the CFHT [5a 
magnitude limits for point sources of u* = 27.4); and • Kg (2.2/im) photometry from the Na- 
tional Optical Astronomy Observatory (NOAO) wide-field IR imager Flamingos (mounted 
on the Kitt Peak 4m telescope) and the Infrared Side Port Imager (ISPI) at the Cerro Tololo 
Inter-American Observatory (CTIO) Blanco 4m telescope. 

The sources in the ground-based catalog were detected from the original (i.e, non- 
matched and non-homogenized PSFs) Subaru i' images, since these had, on average, the 
best seeing condition of ~0'.'6 FWHM, allowing for an optimal deblending of close galaxy 
pairs. Photometric fluxes were instead estimated on PSF-matched images. All magnitudes 
were calculated within circular apertures of 3'.'0 diameter. The final groundbased catalog of 
Capak et al. (2006) includes > 10^ galaxies down to i' = 26 over the entire COSMOS field. 

In deriving the LFs for the COSMOS galaxies, we checked that, in addition to be- 
ing insensitive to slightly different photo- 2; estimates derived by applying or ignoring small 
shifts in the calibration of the groundbased photometry (see Section [5]), the results remain 
unchanged also when using the full ZEBRA Bayesian probability distribution for the photo- 
metric redshift estimates. Furthermore, the LFs do not vary when computing them (a) after 
excluding from the sample all galaxies with ZEBRA values above the 95th percentile of 
the distribution, and (6) by including all galaxies in the sample, independent of the 
value of their ZEBRA fits. 

The ZEBRA photo-z estimates of the COSMOS galaxies are less accurate for Iab > 22.5 
sources not only due to their larger photometric errors, but also due to the lack of photo-z 
calibration at such faint magnitudes. The ongoing spectroscopic follow up of COSMOS on 
the ESO VLT, zC OSMOS, will secure spectroscopic redshifts for > 35'000 COSMOS galaxies 



(ILilly et al.ll2006l ). So far, however, only about 1500 COSMOS galaxies have VLT spectro- 
scopic redshifts, most of which in the 0.2 < z < 1.1 redshift range, at magnitudes brighter 
than lyiB = 22.5. In order to estimate the accuracy of the ZEBRA photometric redshifts 
at magnitudes fainter than Iab = 22.5, we simulated a faint version of the spectroscopic 
redshift catalog by dimming the Iab < 22.5 galaxies with available spectroscopic redshift 
down to our magnitude limit of Iab = 24. We then recomputed the ZEBRA redshifts for 
these artificially-fainted sources. The accuracy of the ZEBRA photo-z's down to Iab = 24 
is a/{l + z) ~0.06. 
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A. 4. Rest-frame quantities 

Total magnitudes for the COSMOS galaxies in our sample were available in the ACS 
/—band. In order to compute the rest-frame absolute -B— band total magnitude (Mb), we 
computed for each COSMOS galaxy the color B — 1^, where B is the observed magnitude 
at Afi^obs = (1 + -2) X ^B, and 1^ indicates the observed J— band magnitude for a galaxy at 
redshift z. The central wavelength of the ACS filter corresponds to the central wavelength 
of the 5— band at redshift z ~ 0.8, therefore B — Iz=o.8 = 0. Defining the color /c— correction 
as K{z) = B — Iz, the rest-frame absolute magnitude Mb can be espressed as: 

MB = h- 5.0 log {dL{z) /lOpc) - 25.0 + K{z) + 2.5 log (l + z), (Al) 

where dL is the luminosity distance at redshift z. We computed the color B—I^ using the best 
fit spectral energy distributions derived for each galaxy as part of the ZEBRA photometric 
redshift calculation. 

In Figure [19] we show 5 — as a function of redshift for our six basic spectral templates 
going from an early-type galaxy template (T = 1, solid line), to the starburst galaxy template 
(T = 6, long-dash dot line). As expected, at redshift z ~ 0.8 the A;— correction is zero for all 
photometric types. 

In the calculation of the LFs, Vmax was computed for each galaxy taking into account 
its ZEBRA photometric type. Specifically, given the galaxy rest-frame Mb absolute mag- 
nitude, we derived the observed expected /—band magnitude as a function of redshift using 
the galaxy best-fit template derived from the ZEBRA photometric redshift. The J— band 
observed magnitude as a function of redshift is required to find ziq and Z24, which define the 
redshift range within which the galaxy could still be seen in the survey, and therefore Knax- 

B. Uncertainties and systematics in the morphological measurements 

To quantify the errors on the measured parameters for the COSMOS galaxies, and to 
assess systematic trends as a function of galaxy brightness, we performed extensive simula- 
tions based on real bright COSMOS galaxies. In particular, we selected 35 bright isolated 
COSMOS galaxies with magnitude between I = 17.5 and / = 18 and, for each of them, we 
created 40 S/N-degraded versions of their images. This was done by dividing the original 
frame by a factor / such that the final magnitude of the object was 2.5 x log(/) fainter 
than the original one. The factor / was chosen so that each galaxy uniformly sampled the 
magnitude range from 18 to 24. The scaled images were then added back to an ACS tile. 
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in randomly-selected positions. Before adding the rescaled image, the original ACS tile was 
multiplied by a factor k, to preserve the noise properties in the final image. Specifically, 
with (Jorig the background noise in the image containing one of the 35 bright galaxies, the 
background noise after the rescaling is cXorig/ /. Adding the rescaled image to the ACS orig- 
inal image increases the noise in the final image to cig^g^i = a^^.^^ + crl^^J P- The factor k is 
therefore computed by requiring that cr\^^^ = a'^^.^^ = crl^^Jk"^ + crl^i^/ p. 

In Figure [2D] we show the results of this set of simulations. In particular we show the 
fractional change of each parameter as a function of the magnitude of the S/N-degraded 
galaxy. Different lines connect the variation of the parameters for each galaxy as a function 
of magnitude. 

For all parameters (except e) the scatter remains low (at a level of 10-20%) and no 
strong systematic effects are observed down to / ~ 22. At magnitudes / > 22 the scatter 
increases and all parameters are systematically underestimated up to a maximum of ~ 20% 
in the faintest magnitude bin. The behavior of e shows a different trend, with the scatter 
starting to increase already at bright magnitudes (at a level of ~ 10%) and systematic effects 
start to be important at / ~ 23. The asymmetry A is systematically underestimated at the 
10% level for /—band magnitudes of order 21, and degrades down to a 20% underestimate 
for Iab = 24. 



B.l. Effects of the Point-Spread-Function 

For very compact objects, the structural parameters computed for the COSMOS galax- 
ies from the HST-ACS images can be affected by the instrumental PSF (FWHM ~ 0.1"). 
In order to determine the limiting size/magnitude at which the effects of the PSF become 
important in our classification algorithm, we performed a set of simulations based on syn- 
thetic galaxy images. Specifically, we constructed artificial images of about 2500 galaxies 
with magnitudes in the range 21 < / < 24, surface brightness profiles described either by a 
deVaucouleur (i.e., = 4) or an exponential {n = 1) profile, and with half-light radii R1/2 
in the range 0.05" < R1/2 < 1.0". These artificial images were convolved with the Tinytim 
(Krist & Hook 1997) HST PSF, and with a PSF constructed from stars present in the COS- 
MOS ACS data. Similar results are obtained with either of these PSFs; in the following, we 
discuss the results obtained with the Tinytim PSF. 

We analyzed and classified, in a similar way as it was done for the COSMOS data, both 
the original and the PSF-convolved artificial galaxies. Comparing the results obtained on 
the original and on the PSF-convolved images we find that galaxies with the exponential 
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[n = 1) surface brightness profile are always classified as Type 2 disk galaxies, and no change 
of class is observed over the entire range of half-light radii and magnitudes explored. For 
n = 4 galaxies, we find that all galaxies are classified as Type 1 (early-type galaxies) when 
the analysis is performed on the not PSF-convolved images. When the convolved images 
are considered, a fraction of about 15% af all galaxies with half-light radius larger than O'.'IT 
changes class from Type 1 to Type 2; 85% of these systems have Iab > 23. For galaxies with 
n = 4 profile, all objects with measured R1/2 < O'.'IT change class from Type 1 to Type 2. 

The fraction of objects in our COSMOS sample with measured R1/2 <0'.'17 is about 4% 
down to our magnitude limit of Iab = 24. For (the steepest-profile galaxies in) this small 
fraction of the sample, the ZEST morphological classification is therefore likely affected by 
convolution with the ACS-HST PSF. 



C. The Kampczyk's SDSS-based z = comparison sample 

The volume of Universe sampled by the Cycle 12 ACS COSMOS data within the < 2; < 
0.1 redshift range is only ~ 5.5 x lO^Mpc^, and includes < 50 galaxies. Therefore, cosmic 
variance and small number statistics in this low-redshift bin are inadequate to provide a 
meaningful 2; ~ calibration for studies of galaxy evolution. 

To study the LF evolution down to 2; = of the ZEST-classified morphological galax- 
ies, we therefore used the SSDS-based dataset of artificial galaxy images constructed and 
presented by Kampczyk et al. (2006). These simulate how the local Universe would be ob- 
served in the COSMOS survey at a redshift z = 0.7. Specifically, the artificial dataset is a 
volume-limited sample of 1813, Mb < —18.5, 0.015 < z < 0.025 galaxies extracted from the 
SDSS-DR4 galaxy catalog. The absolute 5-band agnitude Mb was derived from the SDSS 
g magnitude, by applying a A;— correction computed on the basis of the observed g — r color 



^The SDSS survey provides homogenous color information and spectroscopic redshfits for more than 10^ 
galaxies with magnitude r < 17.77 (Strauss et al. 2002). We expect the Kampczyk et al. SDSS-based sample 
to be unbiased towards environmental effects and fully representative of the local galaxy population down to 
the Mb = —18.5 magnitude limit, since: (a) The photometric+spectroscopic SDSS data are mostly complete 
down to r 17.8 (Strauss et al. 2002), i.e., well below the considered absolute magnitude cut. Specifically, 
the surface brightness cut applied for the selection of the SDSS spectroscopic targets (/i < 23 mag arcsec"^) 
and mechanical constraints of the spectrograph (which exclude galaxies in pairs closer than 55") remove less 
than 1% and '^6% of targets brighter than r = 17.77, respectively. Most importantly, the fraction of excluded 
galaxies does not depend on the galaxy magnitude; and (6) It is extracted from a sample of ~ 700, 000 SDSS 
galaxies over a volume of 2.5 x 10^ Mpc^, which encompasses the full range of local environmental densities, 
from voids to galaxy clusters. 
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and redshift of each galaxy. 

The SDSS (7-band at 2; = almost perfectly matches the ACS F814W (/) band at 
z = 0.7. To simulate how the local galaxies (observed in the SDSS (/-band) would appear, 
in the absence of evolution, in the F814W COSMOS ACS images at a redshift z = 0.7, it is 
thus only necessary to take into account the different pixel scales and PSFs, and cosmological 
surface brightness dimming. The SDSS2=o.7-simulated images do not include either surface 
brightness or size evolution; they were randomly added to the COSMOS ACS images to 
reproduce the ACS-COSMOS noise properties and to simulate galaxy-overlapping due to 
projection effects. 

The simulated SDSS2=o.7 galaxies were analyzed with the same procedure used for the 
real COSMOS data. First, we run SExtractor to detect the sources and to measure their 
/—band magnitudes, position angles and ellipticities. The s ame SExtractor configur ation 



parameters that were used to generate the COSMOS-catalog (ILeauthaud et al. II2006I ) were 
adopted for the extraction of the SDSS2=o.7 objects. About 6% of the SDSS^=o.7 galaxies 
were not recovered by SExtractor; 0.1% were instead detected, but with a magnitude fainter 
than our Iab = 24 limit, and were therefore excluded from the SDSS2=o.7 sample. For each 
simulated galaxy, we then measured the structural parameters described in Section 12.11 and 
applied the ZEST morphological classification described in Section 12.21 

The rest-frame 5— band LF for the constructed SDSS2=o.7 sample is shown in Figure [22] 
[solid histogram). Being based on a volume-limited, random selection of the local galaxy 
population, this LF will have the same shape as the global z = S— band LF. We therefore 
used, to normalize the volume density, the total number of SDSS galaxies with Mb < 
— 18.5, derived integrating the -B— band LF recovered from the (7— band LF of Blanton et 
al. (2003). To convert ^f-band to _B-band magnitudes, we considered two extreme SEDs: a 
young (0.1 Gyr) single burst stellar population with solar metallicity, and an old (12 Gyr) 
single burst stellar population with the same metallicity. The corresponding transformations 
are: B = g + 0.02, and B = g + 0.36. In Figure [22] we show the LFs calculated using the 
normalization derived adopting the "old" SED (dotted line), the "young" SED (dashed line), 
and the average of the two (solid histogram). In discussing the comparison of the COSMOS 
and SDSS2=o.7 LFs, we adopt the SDSS -B— band LF derived with the average SED correction. 
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Table 1. Results of the ZEST PC analysis 





PCi 


PC2 


PCi 


PCi 


PCi 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Eigenvalue 


2.46 


1.19 


0.92 


0.25 


0.17 


Proportion 


0.49 


0.24 


0.18 


0.05 


0.03 


Cumulative 


0.49 


0.73 


0.92 


0.97 


1.00 


Variable 


PCi 


PC2 


PCs 


PC4 


PCs 


Concentration {=xi) 


-0.54 


0.35 


0.18 


-0.34 


-0.66 




0.60 


-0.04 


0.03 


0.39 


-0.70 


Gini {=xs) 


-0.56 


-0.20 


0.14 


0.79 


-0.02 


EUipticity {=X4) 


0.20 


0.57 


0.74 


0.16 


0.26 


Asymmetry {=xs) 


0.02 


-0.71 


0.64 


-0.29 


-0.08 



Note. — Coluiiiiis 2 6 refer to the five principal components 
in ZEST. The first row gives the eigenvalue (i.e., variance) of the 
data along the direction of the corresponding PC. The second and 
third rows show the fraction of variance and the cumulative frac- 
tion of each PC, respectively. The bottom part of the table lists 
the weights assigned to each input variable, in the linear combi- 
nation that gives the direction of the principal component (i.e., 
PCi=Qxi+/3x2+7X3+<5x4+(^X5, with the coefficients o,...,0 given by 
the numbers listed, per each PC, in the bottom part of the table). 
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Table 2. The ZEST classification scheme 





Type 1 
= EARLY TYPES 
(no visible disk) 


Type 2 
= DISK GALAXIES = 


Type 3 
IRREGULAR GALAXIES 


BULGENESS 




from .0 (massive bulge) 
to .3 (bulgeless disk) 




ELONGATION 


from (face on) 
to 3 (edge on) 


from (face on) 
to 3 (edge on) 




IRREGULARITY 


from (regular) 
to 2 (highly irregular) 


from (regular) 
to 2 (highly irregular) 




CLUMPYNESS 


from (smooth) 
to 2 (very clumpy) 


from (smooth) 
to 2 (very clumpy) 


from (smooth) 
to 2 (very clumpy) 


SIZE* 


Rp 


Rp 


Rp 



Note. — Summary of the ZEST classification scheme. Type 1 (early-type galaxies) are spheroids with no visible disk 
(including face-on SO galaxies, for which the identification of the disk component is difficult). Type 2 are disk galaxies, and 
Type 3 are irregular galaxies. A dumpiness parameter is assigned to each (ZEST unit cube and thus) galaxy in the sample. 
Early- type and disk galaxies are assigned an elongation parameter in four steps from (face on) to 3 (edge on), and an 
irregularity parameter (from for regular to 2 for highly irregular galaxies). Type 2 disk galaxies are split in four bins of 
bulgeness parameter, namely 2.0, 2.1, 2.2 and 2.3 from bulge-dominated to bulgeless disks, respectively. Relatively inclined 
SO galaxies occupy cubes classified as T=2.0. 

*The last row indicates that a measure of galaxy size (the Petrosian radius, l lPetrosian Ill976h ') is available for all ZEST— 
classified galaxies, as a byproduct of our analysis. 
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Tablc 3. COSMOS ZEST classification grid 
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Table 3 — Continued 
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Note. — COSMOS-calibrated ZEST grid in each plane of constant value of PC3, 
going from PCz = —2 first panel of the Table, to PCz = 3, last panel. In each panel, 
we give for each PCi, PC2 the ZEST classification as Type 1,2, and 3. For Type 2 
disk galaxies, we also present the "bulgeness" parameter associated with the cube. 
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Table 4. Schechter Function best-fit parameters. 



z(Range) —Ms (Range) Ng^i Mg a 

(Mpc-3 Mag-i) (mag) 
(1) (2) (3) (4) (5) (6) 



All galaxies 



0.2-0.4 


17.0-23.5 


6682 




0.0050 ±0.0006 


-21.03 ±0.25 


-1.26 ± 0.15 


0.6-0.8 


18.5- 


-23.5 


10092 
8146 


T 


0.0049±0.0003 
0.0058±0.0003 
= 1 galaxies 


-21.24 ±0.12 
-21.09 ±0.12 


-1.22 ± 0.10 
-0.95 ± 0.10 


0.6-0.8 


18.5- 


-23.5 


740 
761 


T 


0.00095±0.0001 
0.00095±0.0001 

= 2 galaxies 


-21.07 ±0.20 
-21.12 ±0.20 


-0.26 ± 0.12 
-0.18 ± 0.12 


0.6-0.8 


18.5- 


-23.5 


8820 
7016 


T 


0.0047±0.0003 
0.0053±0.0003 
= 3 galaxies 


-21.04 ±0.12 
-20.89 ±0.12 


-1.25± 0.10 
-0.96± 0.10 


0.6-0.8 


18.5- 


-23.5 


513 
357 




0.0007±0.0001 

0.0005±0.0001 


-20.54 ± 0.4 
-20.27 ± 0.4 


-0.47± 0.20 
0.07± 0.20 










T -- 


= 2.3 galaxies 






0.6-0.8 


18.5- 


-23.5 


3542 
2574 


T-- 


0.0036±0.0002 
0.0031±0.0002 
= 2.2 galaxies 


-20.38 ±0.20 
-20.31 ±0.20 


-1.04± 0.10 
-0.73± 0.10 


0.6-0.8 


18.5- 


-23.5 


3575 
2863 


T -- 


0.0019±0.0002 
0.0024±0.0002 

= 2.1 galaxies 


-20.86 ±0.25 
-20.66 ±0.25 


-1.35± 0.10 
-1.02± 0.10 


0.6-0.8 


18.5- 


-23.5 


1240 
1110 


T -- 


0.0008±0.0001 
0.0009±0.0001 
= 2.0 galaxies 


-21.20 ±0.35 
-21.01 ±0.35 


-1.05± 0.10 
-0.81± 0.10 


0.6-0.8 


18.5- 


-23.5 


463 
466 




0.0007±0.0002 
0.0006±0.0002 


-20.98 ±0.25 
-20.81 ±0.25 


-0.1± 0.10 
0.28± 0.10 



Note. — The columns are: (1) Redshift range: (2) S— band absolute magnitude range; 
(3) Total number of galaxies considered in the fit; (4) Schechter function normalization and 
error; (5) Mg and error; (6) Faint-end slope of the Schechter function. In the z = 0.7 bin, two 
values are provided for each parameter, reflecting the two different estimates for the photometric 
redshifts, derived with (upper values) and without (lower values) corrections to the photometric 
catalogues, respectively (see text). 
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Fig. 1. — Relations between the non-parametric diagnostics (M20, G, A, C, and e = 1 — b/a). 
Contours enclose ~30% (white contour), 80% (grey contour), and 98% (black contour) of 
the galaxies. The main correlations among some of the parameters, such as M20, C, and G, 
are clearly visible in these diagrams. 
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Fig. 2. — Fraction of the total variance explained by each principal component as a function 
of the corresponding principal components for all galaxies down to / = 24 in the ACS-catalog 
(solid circles). Open squares refer to the same analysis performed only on those objects with 
magnitude brighter than I — 22.5. The horizontal hne indicates 20% of the total variance, 
i.e., the value for the five eigenvalues for a sample of 100% uncorrelated variables. 
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Fig. 3. — Illustration of the power of the ZEST grid to separate COSMOS galaxies with 
different structural properties in the three-dimensional PCi, PC2, PC3 space. Shown are 
four separate unit cubes of PC1-PC2-PC3, centered around the values reported in the labels. 
In every unit cube, the few galaxies are representative of the population of objects of that 
specific bin. 
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Fig. 4. — The stamp of a representative galaxy is shown for each cube of PC, arranged 
in planes of constant PC3 in six different panels (a to /). In particular the panels show 
sequentially the PC3 = —2,-1,0,1,2,3 planes. Galaxy structural properties -and thus 
galaxy types- change smoothly through the PC1-PC2-PC3 space. 
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Fig. 4 — Continued. 
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Fig. 4 — Continued. 
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Fig. 5. — Schematic view of the ZEST morphological classification. Shown are slices of 
the PC space at constant value of PC?,. Each slice shows PCi vs PC2; the value of PC3 
increases from bottom-left to upper-right. In each unit cube the symbol indicates the ZEST 
classification associated with that cube (according to the legend shown in the bottom-right 
corner of the Figure). In particular, stars represent T = 3 galaxies, and single ellipses 
represent T — 1 galaxies. The double eUipses indicate T — 2 galaxies, for which the size of 
the inner ellipse increases for increasing bulgeness parameter. 
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Fig. 6. — The fraction of objects with a given RC3 classification (E, SO-Sab, Sb-Scd, and Sd 
and later) that have respectively T=l, 2.0,2.1,2.2,2.3, and T=3 ZEST morphological Type. 
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Fig. 8.— The solid histograms show the distribution of (A(PC)) = Y.i \PCi - PCi,o|)/3, 
with i = 1, 2, 3, PCi^Q the initial values of the PCs, and PCi the measured PCs after signal- 
to-noise image degradation of a set of 17.5 < / < 18 COSMOS galaxies. We show four 
different magnitude bins for the S/N degraded galaxies, from / = 22 (top- left panel) down 
to / = 24 (bottom-right panel. The hatched histograms show the (A(PC)) distribution for 
the galaxies with AT ^ 0. 
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Fig. 9. — Three examples of how galaxies move in the PC space as the signal-to-noise of the 
galaxy image decreases. The start- and end-coordinates in PC space are indicated (PCin, 
and PCfin, respectively). 




Fig. 10.— Distribution of our ~56000 Iab < 24 COSMOS galaxies in the G-M20 plane (top 
panels) , C — M20 plane (middle panels) , and C — A plane (bottom panels) . In each column 
we highlight the position of a different ZEST morphological class by drawing the contours 
enclosing 99% of the objects in that ZEST class. In the first, second and third columns we 
show respectively T = 1, 2 and 3 galaxies (red, blue and green contours). The blue color for 
T = 2 galaxies ranges from dark blue for bulge-dominated T = 2.0 galaxies to light blue for 
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Fig. 11. — /—band magnitude distribution for the ACS-galaxies with no detection in the 
ground-based catalog (top panel), and for the ACS-galaxies whose ground-based match was 
associated with more than one ACS-detected galaxy. Both distributions are normalized to 
the magnitude distribution of all galaxies in the sample. The fraction of missed objects stays 
relatively constant down to Iab ~ 22 and then increases toward fainter magnitudes. 



-53- 



I 



to 
I 



-ad 
-5J 



-ad 



-4J 



! f ALL' 




























1 ■ 





— iifc— stLfc- _iit_a(Lfc_ifii _iit_a(Lfc_ifii _iit_a(Lfc_ifii 



M 



B 



Fig. 12. — The 1/Knax corrected (stars) and uncorrected (dotted line) LFs derived using the 
ZEBRA Maximum LikcUhood photometric rcdshifts for the global sample (top row), for T =1 
early-type galaxies (second row), for T =2 disk galaxies (third row), and for T =3 irregular 
galaxies (bottom row). In each row, we show the LFs in four rcdshift bins: 0.2 < z < 0.4 (first 
column), 0.4 < z < 0.6 (second column), 0.6 < z < 0.8 (third column), and 0.8 < 2; < 1.0 
(last column). Error bars in each luminosity bin take into account Poissonian errors only. 
In each redshift bin LFs are plotted down to the S— band at which the sample is complete, 
regardless of the color of the galaxy. Sohd curves show the best-fit Schechter function to the 
COSMOS LFs in the redshift interval 0.6 < 2; < 0.8, (de-)brightened in each redshfit range by 
~ 1.3 X z magnitudes (see text). The stars show the LFs obtained using the ZEBRA photo- 
2;'s derived after correcting the photometric catalogue for detected systematic calibration 
errors; grey shaded regions show the uncertainty in the LFs due to the use of photo-^'s 
derived without applying any correction to the photometric catalogue. 
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Fig. 13. — Results of the simulations performed to assess the impact of the photometric 
redshift errors in the calculation of the LF, in four redshift bins, as indicated in each panel. 
Open diamonds represent the median volume density computed from the 100 mock realiza- 
tions; the associated grey area indicates the 16*'*-to 84*'* percentiles of the distribution in 
each magnitude bin. Points are shifted of 0.05 magnitude from the center of the magnitude 
bin for presentation purposes. The solid line with errorbars represents the LF of the original 
mock catalog used to generate the 100 mock datascts. For comparison, we also show the LFs 
derived using two independent mock catalogs (dashed lines). These LFs show that cosmic 
variance in the lowest redshift bin dominates the uncertainty budget in the calculation of 
the LF errors. 
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Fig. 14. — Corrected 1/Knax LF computed for for T = 2 galaxies, split by their bulgeness 
parameter, from T = 2.3 bulgeless disks (top row), to T = 2.0 bulge-dominated galaxies 
(bottom row). Dot-dashed curves in the redshift bin 0.6 < z < 0.8 show the LFs derived for 
the SDSS^=o.7 comparison sample. The solid curves show the best-fit Schechter function to 
the COSMOS LFs in the redshift interval 0.6 < z < 0.8, brightened in each redshfit range 
by ~ 1.3 X 2 magnitudes (see text). Grey shaded regions are as in Figure | 
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Fig. 15. — Comparison between the corrected COSMOS LFs (stars) and the (redshifted) 
SDSS2=o.7 LPs, representing the no-evolution predictions at redshift z = 0.7 (dash-dotted 
hne). The sohd curves shows the SDSSj;=o.7 LFs brightened by 0.95 magnitudes (see text for 
details) . Grey shaded regions are as in Figure [121 
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Fig. 16. — The total radius Rtot versus /—band magnitude for all COSMOS galaxies with 
Iab < 24. All objects with CLASS_STAR> 0.6 are identified as solid circles. Down to 
Iab = 24, stars form a tight sequence in this plane, which can be identified using the 
SExtractor CLASS_STAR parameter. 
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Fig. 17. — Comparison between the Petrosian radii of Iab < 24 COSMOS galaxies, as 
derived in our analysis, and the half-hght radii derived by Sargent et ah (2006) for the 
Iab < 22.5 subsample (using GIM2D fits to the galaxy surface brightness distribution). The 
solid and dashed lines represent the relation between Rp and R1/2 expected for a Sersic 
profile with n = 1 and n = 2, respectively. 
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Fig. 18. — Examples of the steps performed to compute the Petrosian radius Rp. In particu- 
lar, we show four galaxies with increasing magnitudes, from I — 18.80 (top row) to / = 22.3 
(bottom row). For each galaxy we show the original stamp (first column), the cleaned stamp 
(second column), and the segmentation map (third column). In the last column we show 
the surface brightness profile (top panel) and the ratio rj = fi{R) / {fi{< R)) (bottom panel) 
as a function of R in arcseconds (solid circles). The solid horizontal line in the rj — R plots 
shows r) = 0.2. In the last example, rj decreases to a minimum value larger than r] = 0.2 and 
then increases again. This behavior is due to contamination to the external isophotcs from 
light coming from the wings of the nearby source, still visible in the cleaned stamp, despite 
the attemp to remove it before proceeding with the analysis. This hght causes the surface 
brightness profile to fiatten at large radii, as seen in the ji — R diagram. To correct for this 
effect, we subtract the residual light contribution from the nearby source. The correction is 
derived by fitting a linear relation (shown as dashed line in the — i?plot) to the external 
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Fig. 19. — Color /c— correction (defined as k{z) — B — I^, see text for details) as a function of 
redshift for different spectral templates, going from an early-type galaxy template (T = 1, 
solid line), to the starburst galaxy template (T = 6, long-dash dot line). At redshift z ~ 
0.8 the color A;— correction is virtually zero for all photometric types, since at this redshift 
the rest -frame 5— band coincides with the /—band filter from which the 5— band absolute 
magnitudes are derived. 
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Fig. 20. — Results of the simulations performed by degrading the S/N of real galaxies (see 
text for details). Lines in each panel show the fractional variation of the parameter as a 
function of magnitude (each line shows a different galaxy). From top to bottom we show 
ellipticity, concentration, Gini, M20, and asymmetry. Systematic effects are < 20% over the 
whole range of observed magnitudes. 
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Fig. 22. — The rest-frame 5— band LF for the SDSS galaxy sample of Kampczyk et al., which 
we use to quantify the evolutionary trends since z 1 pres e nt in our COSMOS sample (solid 



histogram). This is normalized using the iBlanton et al. I (120031) (7— band LF, converted to 
the 5-band magnitude using the average transformation (solid line) between a 12Gyr-old 
single burst SED with a solar metallicity (dotted line) and a single-burst O.lGyr-old SED 
with a similar metallicity (dashed line). 
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